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Abstract 

We relate Nekrasov partition functions, with arbitrary values of ei,e2 parameters, to 
matrix models for /3-ensembles. We find matrix models encoding the instanton part of 
Nekrasov partition functions, whose measure, to the leading order in €2 expansion, is given 
by the Vandermonde determinant to the power /3 = —ei/e2- An additional, trigonometric 
deformation of the measure arises in five- dimensional theories. Matrix model potentials, 
to the leading order in £2 expansion, are the same as in the /3 = 1 case considered in 
0810.4944 [hep-th]. We point out that potentials for massive hypermultiplets include multi- 
log, Penner-like terms. Inclusion of Chern-Simons terms in five-dimensional theories leads 
to multi-matrix models. The role of these matrix models in the context of the AGT 
conjecture is discussed. 
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1 Introduction and summary 



Quantum gauge theories with M = 2 supersymmetry appear to possess a rare feature 
of getting more and more fascinating as the time flows by. Since the seminal solution of 
Seiberg and Witten jll [2], it is known that their solution is encoded in Seiberg-Witten 
curves. The very existence of these curves inspired several important developments and 
led to remarkable relations to other branches of physics and mathematics. 

The first such branch is a theory of random partitions. In [21 HI O |6l [7] it was shown 
that Seiberg-Witten curves arise as limiting shapes of large ensembles of random two- 
dimensional partitions. These partitions label instanton configurations, which are the only 
contributions to the gauge theory partition function computed by localization. To apply 
localization techniques one has in fact to consider a more general setup, which, among 
several other subtleties, involves a non-trivial spacetime metric called the f2-background. 
This background is characterized by two parameters ti and €2- Ultimately one can explicitly 
write down partition functions of gauge theories Z, which now bear the name of Nekrasov, 
and in particular they depend on ei and €2- The information about the Seiberg-Witten 
curve is encoded in the prepotential Fq, which is the leading term of the free energy 
Fq = — lim^^ ,=2^0 £1^2 log ^- In the special case ei = —€2 = h the free energy has the 
expansion 

00 

F = logZ = 5^F,^2^-2, (1) 

9=0 

and higher genus terms Fg are also associated to the ordinary Seiberg-Witten curve. For 
arbitrary €1,62 this expansion takes more general form and it is natural to expect that 
it should be associated to some deformation of the Seiberg-Witten curve. As we explain 
further there are strong indications that the non-commutative deformation arises in such 
a general case. 

Another class of systems, whose solutions are encoded in complex curves, are matrix 
models, defined by matrix integrals 

Z = J VMe--^^'''^^'\ (2) 

The integration measure can typically be diagonalized and expressed in terms of eigen- 
values Xi of matrices M. In case of hermitian matrices, such diagonalization results in 
the Vandermonde determinant VM = Yli dx-i Y\i<,j{^i ~ ^jY- The solution of such matrix 
integrals can also be expressed in the form ([1]), which is encoded in the underlying spectral 
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curve. One is therefore tempted to devise such (hermitian) matrix models, whose spectral 
curves would be identified with Seiberg-Witten curves, and moreover the full solutions 
([1]) on both sides would match. Such a program was initiated by Dijkgraaf and Vafa in 
[SI |9], where they related both systems by a chain of dualities which involved topological 
string theory. The connection to topological strings is not a surprise: their amplitudes 
also take form ([T]), and it is known that M = 2 theories can be engineered by considering 
topological strings on non-compact Calabi-Yau manifolds [10], whose geometry is encoded 
in yet another complex (B-model) curve [11]. This relation is exact, i.e. topological string 
partition functions agree with Nekrasov partition functions to all orders in genus expansion 
[121 [13] • This agreement was also extended, in terms of refined topological vertices [HI US], 
to arbitrary values of ei, €2 in |T6| [T7] . 

However, the relation of Seiberg-Witten (and therefore topological strings) to matrix 
models is in fact more subtle: while it was shown that the expansions ([T]) agree for the 
lowest order terms -Fo)-^i even for simple polynomial potentials [Ml IS], the agreement of 
the full series ([T]) is much harder to achieve and depends on details of matrix potentials 
V{M) [20]. Nonetheless, recently matrix models have been found [2T| |22] . which by the 
very construction agree with the Nekrasov partition functions, and therefore automatically 
lead to the identification of full expansions (pQ) on both sides. 

Recently the above web of dualities has been extended again, in a novel and surprising 
direction: it was observed that Nekrasov partition functions agree with correlators in 
Liouville theory [23j. In view of the connections described above, one might therefore 
expect that there are connections between Liouville theory and matrix models. Such 
connections have been proposed recently by Dijkgraaf and Vafa [21], in a way which extends 
their earlier ideas mentioned above. We will discuss these relations more in what follows. 

In the above web of dualities matrix model possess one peculiar feature. While the 
amplitudes in gauge theories, topological strings]^ and Liouville theory depend in a natural 
way on the two parameters 61,62, it is not a priori clear how (or whether at all) the 
dependence of the matrix integral ([2]) on /I can be refined. Nonetheless, matrix models 
are essential ingredients of all those dualities, and therefore they would appear incomplete 
without such a refinement. 



^Topological strings, as such, are defined for a single parameters gs = (1 = —£2- However their 
amplitudes can be refined in a natural, combinatorial way to two parameters using the construction of 
refined vertices |14[ I15|. It is argued in [T7j, that for topological strings on toric Calabi-Yau manifolds 
which engineer gauge theories, the amplitudes computed using these two vertices give the same results. 



4 



In this paper we show that such a refinement of matrix models indeed exists, and 
it takes form of the so-called /3-deformation. By explicit computation, we transform a 
large class of Nekrasov partition functions into the form of /3-deformed matrix models, i.e. 
matrix models of the form ([2]), however with modified determinant [25] [261 127]. the four 
dimensional case this determinant is given by 



and the exponent /3 captures the dependence on parameters of the ^-background 



Clearly the special case ei = —62 = ^, or /5 = 1, corresponds to the ordinary Vandermonde 
determinant. We derive /3-deformed matrix models for four- and five-dimensional U{n) 
theories with massive hypermultiplets, as well as five- dimensional Chern-Simons terms. In 
the five- dimensional case the determinant ([3]) is in addition deformed in the trigonometric 
way. The generalization to six dimensions is straightforward, and results in the elliptic 
deformation, similar as in [13| 122] . The methods which we use are essentially extension of 
those introduced in [2T| [22] |28] : for related developments see also |29[ \3U\ l3T| [32] . 

It was already argued in |24) that refinement of matrix models should take such a 
form, based on properties of the Q background, as well as the example of the /3-deformed 
Gaussian matrix model. This Gaussian model is a special case of the Selberg integral, and 
its exact value is known. However this example seems to relate only to four-dimensional 
theories, and it is not clear how to extend it beyond the quadratic potential. 

More precisely, the /3-deformed matrix models in the form given above, arise in our 
computation in the leading order in the expansion in one of the parameters, say €2- The 
/3-deformation of the measure is extracted from the asymptotics of its yet more general 
form. In addition there is the whole sub leading series in 62, both in the measure, as well as 
in the matrix model potential. This is analogous to the results in [211 [221 [2S1 [21] for /3 = 1, 
where matrix potentials (but not the measure) included additional subleading series in h. 
An additional subtlety in the present /3 7^ 1 case is that these subleading terms cannot 
be immediately symmetrized (in a sense which will be explained in what follows), which 
makes their matrix model interpretation less clear. Nonetheless, we postulate that these 
subleading terms are essentially inessential, and the knowledge of just the leading terms 
should be sufficient to analyze the entire theory. This is so in matrix models corresponding 
to /3 = 1 analyzed in |21[ [221 [28] . A detailed explanation of this phenomenon is given in 




(3) 




(4) 
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|29| . where analogous matrix models for Hurwitz numbers are derived. In these cases the 
leading piece of the matrix model is sufficient to solve the whole matrix model, because it 
encodes the spectral curve, and the knowledge of the spectral curve is sufficient to define 
recursively symplectic invariants of Eynard and Orantin [33j . In case of curves arising in 
the B-model topological strings, these symplectic invariants agree with Gromov-Witten 
invariants [Ml ES]- In fact, the results of ItSSj have been recently extended to the matrix 
models for /5-ensembles [261 [27]. We therefore suppose that similar scenario should hold in 
the present case. 

The organization of this paper is as follows. In the reminder of this section we sum- 
marize in detail our results, and discuss their relation to the AGT conjecture. In section 
|2] we present Nekrasov partition functions in the form appropriate for our purposes. In 
section [3] we present a general scheme of deriving matrix models associated to Nekrasov 
partition functions. In section H] we provide such a derivation for four- dimensional the- 
ories. In section |5] we provide a derivation for five-dimensional theories. In section [6] 
we discuss directions for further studies. In appendices we present yet more details on 
Nekrasov partition functions, explain manipulations with infinite sums, present various 
asymptotics necessary to extract /3-deformed measures, and provide a complete example 
of the subleading terms in the four-dimensional matrix model potential. 

1.1 Summary of the results 

In this paper we find /3-deformed matrix models, which encode Nekrasov partition 
functions for theories U{n) gauge groups, with various matter contents and in various 
dimensions. We find it convenient to rescale the constant in front of the potential in ([2]), 
and write the /3-deformed matrix models in the form 

Z= [ VMe-^^'''^^\ (5) 

where M G Mat^N denotes nN x nN matrices from /3-ensembles. 

For four-dimensional theories, to the leading order in €2, we find that the measure DM 
involves the /3-deformed Vandermonde determinant 

VM = \{dx,\{{x,-Xjf^ (6) 

i i<j 

while four-dimensional potentials take form 

V'\X) =tX + KtKx) + V^lu)funAx)- 
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The linear term tx encodes a dependence on the instanton counting parameter, while the 
contributions from vector and hypermultiplets in the fundamental representations, to the 
leading order in €2, are given by 

n 

KVc{x) = 2 ((a; - ai) log(a; - ai) - {x - a;)) , 
1=1 

Vfundi^) = [ - (x - TT^f) ^og{x - mf ) + (x - mf . 
f=l 

The expression for V"^tfti/Mnd(^) analogous to V^^^^{x), but with masses mf replaced 
by ei + €2 — rrif. We note in particular that potentials for matter in ( ant i) fundamental 
representation include Penner-like factors of the form 

rrif log{x — rrif). 

f 

Analogous factors were obtained from completely different perspective by Dijkgraaf and 
Vafa in [21]. While the arguments of logarithms in were chosen, at least to some 
extent, at will due to conformal invariance, in our case they are fixed as (x — rrif). 

We also derive matrix models for 5-dimensional theories compactified on a circle of 
radius R^. We usually suppress the dependence on this radius, and it can be reintroduced 
by rescaling x and other parameters (such as a;,mf) by R^. The matrix models for five- 
dimensional theories also take form (|5]), with the measure which involves, to the leading 
order in €2, the Vandermonde determinant with both trigonometric and /3-deformation 

VM = l[ dx, n (2 sinh (7) 

i i<j 

The five-dimensional potentials read 

V'ix) =tX + Vliix) + V^liu)funA^)- 

The linear factor tx encodes in particular the dependence on the instanton counting pa- 
rameter, while contributions from vector and hypermultiplets, to the leading order in 62, 
are given respectively by 

n 

KM{x) = ^x^ + 2 5^Li2(e--+»0, 

1=1 

V^U^) = 5^Li2(e^-™0- (8) 
f 
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The contribution for antifundamental matter V^^^^j^^j^{x) is related to Vj^^^{x) as in the 
four-dimensional case. 

In five dimensions one can also introduce Chern-Simons terms, which are parametrized 
by a single integer mcs- In this case we obtain multi-matrix models with additional linear 
terms 

n N 
1=1 i=l 

where xf^ represent eigenvalues of the Vth set of matrices. 

We stress that all these potentials are the same as in the /3 = 1 case discussed in pT] . 
In the present case the only modification arises as the /3-deformation of the Vandermonde 
determinant. 



1.2 Relation to the AGT conjecture 

As we already mentioned, in [23] remarkable connections between four-dimensional 
Seiberg-Witten theories and two-dimensional conformal field theories have been conjec- 
tured. This observation was again motivated by properties of Seiberg-Witten curves: in 
[36j they were related to other curves, whose gluing allows to construct more complicated 
gauge theories from simpler ones, in a way analogous to the construction of conformal field 
theories on arbitrary Riemann surfaces. Then direct comparison of certain quantities in 
several theories from both these classes led to general conjectures. These quantities involve 
one- loop and instanton parts of Nekrasov partition functions Z^~^°°^, ^I'^nst ^y^q Seiberg- 
Witten side, and Liouville three-point functions C^^^^ given by the DOZZ formula and 
conformal blocks J-'fl on the CFT side. 

In particular, under this correspondence, the central charge of the Liouville theory 

62 

is expressed in terms of the ratio of ei and €2- In the special case ei = —€2 = h we obtain 
c = 1 and the Liouville theory reduces to the free fermion theory. Further aspects and 
checks of this correspondence were discussed, among the others, in [371 EHl EH SOl HH I32| 

sal El sa SSI mi SUSS]. 

Dijkgraaf and Vafa proposed a derivation of this conjecture inspired by topological 
string dualities and relation to matrix models [23]. A sequence of these dualities is shown 
in figure [U The top arrow represents the original AGT conjecture |i23j. The bottom arrow 
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Figure 1: Dualities relevant to the AGT conjecture. The top grey arrow repre- 
sents the original AGT conjecture f^ . The other grey arrows represent a chain 
of dualities proposed by Dijkgraaf and Vafa ^24j- The red arrow represents our 
present explicit reformulation. 



refers to general relations between topological strings and matrix models discovered by 
Dijkgraaf and Vafa in |8l|9]. The crucial ingredient of their recent derivation [21] involves 
engineering Calabi-Yau manifolds relevant for theories considered in [22] , and is represented 
by the vertical arrow on the left. The chain of dualities is completed by the arrow on the 
right, representing the relation between matrix models and conformal field theories [SOtlST]. 
In fact, to make this chain of dualities rigorous, several issues are still to be completed. 
For example, the relation between matrix models and conformal field theories is known to 
hold rigorously only in c = 1, or equivalently (3 = 1 case. Also the relations to topological 
strings discussed in [24J were provided for /3 = 1 case, and their extension beyond this case 
is more subtle. Let us also note, that various aspects of the matrix models proposed in 
[21] were further analyzed in 03 SHI 113 SHI ©J • 

We wish to point out, that our present results can be regarded as an explicit realization 
of a part of the program proposed in [2^. Our results are represented by the red arrow in 
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figure [U we directly relate Nekrasov partition functions to /3-deformed matrix models. As 
argued earlier, we suppose that the leading form of these matrix models should be sufficient 
to solve entire theory, due to powerful recursive relations found in [26 | l27 t 133 ] . (In fact, due 
to equivalence of five-dimensional gauge theories and topological strings |12[ [T3| [T6| [T7] . 
our results automatically realize the equivalence represented by the bottom arrow, which 
therefore could also be drawn in red.) 

Furthermore, the relation to Liouville theory might be completed due to the equivalence 
of non-commutative (or quantum) curves arising on both sides. It was proposed in [23] how 
to associate a quantum curve to the Liouville theory. This curve is given by the operator 
equation 

x' = T{z), 

where T{z) is the Liouville energy-momentum tensor, and it is supposed to encode the data 
of the Liouville theory. On the other hand, the non-commutative spectral curves have been 
recently associated to /3-deformed matrix models in [261 [27], ^ind it was shown how to gen- 
eralize to this case the recursion relations familiar from the (3 = 1 case [33]- Even though 
matrix models considered in [26l [27] have only polynomial potentials, one might hope to 
extend those results to the potentials which we obtain here. In such case, analogous non- 
commutative curves might be associated to our matrix models. The adventage of such an 
approach would be twofold. Firstly, the corresponding non-commutative symplectic invari- 
ants should hopefully agree with refined Seiberg-Witten and (for five-dimensional theories) 
refined Gromov-Witten invariants. This would provide more rigorous proof of the AGT 
conjecture. Secondly, this would provide a concrete proposal for the non- commutative gen- 
eralization of the Seiberg-Witten curve - i.e. if such a generalization exists, it is natural 
to expect that it should agree with the (non-commutative) matrix model curve, similarly 
as in /3 = 1 case |^34J- It would be very interesting to understand such a generalization 
of the Seiberg-Witten curve purly from the gauge theory perspective, as well as from its 
embedding in string theory. In this regard we note that yet another form of quantum 
curves have been considered in the context of Seiberg-Witten theories in [52| [53] . which 
have a string theory interpretation in terms of intersecting D4- and D6-branes. Suppos- 
edly those curves should also be related to the non-commutative curves mentioned above, 
and the B-field present in the intersecting brane system would provide physical reason 
for the non-commutativity. The relations between matrix models and non-commutative 
curves are discussed also in [45j. Connections between all these issues are currently under 
investigation. 
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2 J\f = 2 theories and Nekrasov partition functions 



Partition functions for M = 2 theories were derived by Nekrasov using localization 
techniques in [3] and analyzed further in great detail, among the others, in [U [5l El [TTJ 
[TB| [25] . They consists of the one-loop and instanton part. In this paper we consider only 
instanton parts, and refer to them simply as (Nekrasov) partition functions. 

Partition functions for U {n) theories can be viewed as a generalization of the Plancherel 
measure to the case of n sets of partitions j^l EI]. They are given by a sum over a set of 
n partitions Xf\ with I = 1, . . . ,n labeling various partitions and index i > denoting 
the i'th row of a given partition. These partition functions depend the scalar vevs ai via 
o;fc = CLi — dk, the dynamically generated scale or other counting parameter A, and the 
parameters of the f2-background ei,e2- We also introduce 



which will play an important role from the matrix model perspective. Of course a de- 
pendence on masses mf arises in theories with A^f massive hypermultiplets. We introduce 
also 



In five-dimensional theories there is also a dependence on the radius of the fifth dimension 
i?5, which we usually suppress. It can be reintroduced by rescaling quantities which appear 
in our final expressions by R^. 

The tuples of partitions label various instanton configuration, and from the localization 
perspective the corresponding contributions are naturally expressed in terms of so-called 
arm-lengths and leg-lengths of boxes in these partitions |3l HI [5l [6]. These expressions are 
summarized in the appendix [XI However these expressions can also be written purely in 
terms of the length of rows of these partitions, which is more useful from the perspective 
of putative matrix models we are after. As reviewed in the appendix 13 in this form the 
instanton partition functions for vector and (anti)fundamental multiplets read 




(9) 





Z 





A=(A{i),...,A(")) 



Zi' 




1 





(10) 



2n|A| 



T{X? - + /3(j - z) + kk) r(/3(j + kk + f3) 
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where |A| = Yll^=i X]j l-^i'^l- The contribution to ^I'lj^^jj^^^ is analogous as to Z^^^^^, but 
with TJif replaced by ei + 62 — mf . 

In the five- dimensional case we introduce 

q = e'\ t = e-'\ Qi = e""' = e'^^\ Qi^^ = e'''""' = e'^^'K (12) 

The partition functions read (again see the appendix 

/ J \,vec \,{anti)fund \;tncs'' 

A=(A{i),...,A(")) 

erl\,ii,)(Q,,g^-'''-ft.-^-i;g)^ (Q,.t.-;g)^ 

, = TTTTTT^I (14) 

where ||A|| = Yli-^f- The contributions for hypermultiplets in antifundamental represen- 
tations Z'i'^ arise analogously as in the four-dimensional case. The term Z^'^ 

\,antifund \,mcs 

encodes Chern-Simons terms which may arise in the five-dimensional case, and which are 
parametrized by a single integer mcs'i such Chern-Simons contribution vanish in the four- 
dimensional limit. 

More generally, one can also consider contributions from matter in bifundamental or 
adjoint representations. The generalization of matrix models presented here to those cases 
is not difficult, and we leave it for future work. In principle, extending our methods to 
bifundamental matter would lead to multi-matrix models, while analysis of adjoint matter 
would generalize the results of [22j. 

We therefore focus on contributions for vector and (anti)fundamental multiplets, as 
well as Chern-Simons terms. To sum up, we note that the partition functions given above 
can be written, up to inessential constant factors, in a unified way as 

7 — A2n|A|7 y yCS 

^ / J ^vec^{anti)fund^rn,csi 

A=(A{i),...,A{")) 

^ ^ 1-r /(aP - Af ) + /3(j - ^) + + f5) f{k, + /3(j - ^)) 

- mm /(6,-M,-./3+i) ('^^ 
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where, respectively in the four- and five- dimensional case, the function f{x) takes the form 

f'^ix) = r(x), fix) = , \ ■ (18) 

In five dimensions the Chern-Simons term Z^^^ can also be expressed purely in terms of 
the lengths X^p, which we will discuss in more detail in section [5l In the four-dimensional 
case we simply put Z^^^ = 1. 



3 Matrix model representation 

Nekrasov partition functions involve sums over n-tuples of partitions A^'-*. In the pre- 
vious section we expressed these sums purely in terms of the lengths of rows of these 
partitions Xf\ In this section we reinterpret these lengths as eigenvalues of matrices in 
certain ensembles. We follow the general strategy presented in [2T | |22 | [28]. In this section 
we do not consider Chern-Simons terms in five- dimensional case, or assume that rrics = 
in (IT^ . Arbitrary mcs will be considered in section [521 

To start with, we truncate the sums over partitions in Nekrasov partition functions, to 
sums over partitions with at most non-zero rows. After several manipulations we will 
find a matrix model representation of Nekrasov partition functions in the large A^ limit of 
ensuing expressions. Having fixed A^, we introduce 

h? = -t + N + k, (19) 

as shown in figure |2J In terms of these variables, the expressions (I16p and (II 7p take 
respectively the form 

^ ^ i-r f{h? - hf + j)(l + f{h, + /3(j - .)) 



- linn jib^rj^;^wTT) • ^''^ 

The first important observation is that the above expressions can be written in terms 
of a single set of variables obtained from concatenation of sequences 

h=l,...,nN = {h^i \ ■ ■ ■ , h^-N 1 ' ^i^'*' • • • ' ^7V^)- (22) 

Assuming that a/'s are decreasing and in sufficiently large distances from each other, the 
sequence li=i^,„^nN is decreasing. As shown in appendix [B| the expressions (120|) and (12T|) . 



13 



,m 1.0) An) ,{n) ,(«)_, 



Figure 2: Translation of 's to h^P 's and li 's, which will be reinterpreted as 
matrix eigenvalues. The leg-length and arm-length of given box (in black) are 
defined respectively as the number of boxes to the right, and above this box (in 
yellow). 



for arbitrary / (not necessarily of the form f|T8|) ). can now be written as 



nN 



^vec 



( f{(3) \ fr f{h - I, + i^modN - JmodN)il ' + ^ 

i^^ f{-{h-k + {ZmodN-N){l-P))) 

^ LLLL _h. ^{^ ... _ AT\(T _ R\ M 

1=1 i-- 

Nf nN 



}^ f{k-bl + itmodN-N)il-P)+(3) 

- Mf- N + imodN{l- P, 
fiPi - Mf - IrnodNP + 1) 



^fund - 1111 771— 17— — --^Tn ^^^^ 



f=l i=l 

where imodN = imodN. Moreover, as k are expected to be large, we introduce rescaled 
variables 

Xi = e2li. (26) 

Now we wish to reinterpret the above expressions as partition functions of /3-deformed 
matrix models, with Zj's playing the role of eigenvalues. As we show below, the functions 
f{x) indeed have a form relevant for such an interpretation, at least in the leading expansion 
in £2- As we will see, the subleading terms in such an expansion depend not only on 
variables /j, but also imodN . This makes matrix interpretation of those subleading terms 
less clear, however, as explained before, they should be inessential to the solution of our 
matrix models. We also note, that such dependence on i in the subleading terms does not 
arise in the /3 = 1 case, see also |2H 122] . 
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To relate the above expressions to matrix models, in what follows we will perform the 
following steps: 



1. Reinterpret the factor (125]) in Z^_^^^ involving differences (/j — Ij) = (xj — Xj)/e2 as 
the (possibly deformed) Vandermonde determinant 

2. Reinterpret other factors (I24p . (125 p . etc., involving just a single li = Xi/e2, as contri- 
butions to the potentials Kec? Vfund, etc., by writing them in the form 



nN 

r i 

exp 



^ nN 



£2 . , 

1=1 

3. Introduce an auxiliary function fpoies{x), which has simple poles at all integer values 
of the argument 

4. Replace the summation over sequences of li by the sum over all k using symmetriza- 
tion, and subsequently by the integration over Xi = €2(1 over a contour C which 
encircles a (part of) real axis: 

nN 

{Xi 



(28) 



Then, in the leading expansion in 62, the resulting expressions have the form of the eigen- 
value representation of matrix integrals. In the following sections we discuss this program 
separately for four- and five-dimensional theories. 



4 Matrix models for four-dimensional theories 

In this section we reinterpret four-dimensional Nekrasov partition functions as matrix 
models in the leading expansion in €2, and discuss the example of SU{2) theory with two 
fundamentals and two antifundamentals. 

4.1 Derivation of matrix models 

Here we follow the steps 1 — 4 listed in the previous section, and apply them to four- 
dimensional Nekrasov partition functions. 
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1. The factor (123 p for the four-dimensional theory is reahzed in terms of f^'^{x) = T{x) 
functions (ITS]) . This factor is therefore a ratio of gamma functions. Here comes the 
crucial observation: the relation to the /3-Vandermonde arises from the asymptotics 
of the ratio of gamma functions (156|) presented in appendix IC.2I As we expect 
to be large, we can make the following identification between and the arguments 
of f^^ 

Z = li- Ij, a = {imodN - jmodN){l " /?), /? = (29) 

Therefore, expressing this expansion in terms Xi introduced in (l26l) . we get 



TT r(/i - Ij + {imodN - jmod Af)(l - (3) + (3) ^ ^_/3„iv TTT/ _ ^ ( I -I- Qi — 

^(/,-/, + (^^„,;v-J„^od7v)(l-/3)) l}^""' V ^ 

As presented in appendix IC.21 the subleading terms are of the form 

OO 1. 

o.)-, rr- 



(30) 



nr. - — T • ^ 



We note that the dependence on the coefficient a = {imod n — jmod n){^~ P) is encoded 
only in higher order coefficients C„(/3,a). Note that these coefficients vanish in the 
(3 = 1 case. This is why no ^-deformation of the measure was observed in |2H |22] . 

2. Using the expansion of the logarithm of the gamma function fl53p . the remaining fac- 
tors from Z^^^, given in (12^ with f^'^{x) = T{x), can be interpreted as contributions 
to the potential in (127|) . The full form of Vjf^{x) is given in (158|) . To the leading 
order in €2 it reads 



Ktc(^) = 2 V (x - aO log(x - ai) -{x- ai) + 0{e2) . (31) 



n 

;=i 

The same potentials were obtained in the leading order in the /3 = 1 case in |21] . As 
explained earlier, such leading contributions should be sufficient to get the spectral 
curve of the matrix model, and subsequently solve it. In fact the subleading terms, 
for general (3, do not depend only on A^^ but also explicitly on (imodN), therefore 
they cannot be simply symmetrized and their matrix model interpretation is less 
clear. Nonetheless, as an example we present the full form of V^^^{x) in appendix [Dl 
Similarly as above, we also get the contribution for fundamental multiplets 

Nf 

V}Ld{x\ m) = ^ [ _ (x - mf ) log(x - mf ) + {x- m,) + ^(es)] , (32) 

/=i 
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and we omit the constant factors of the form Ylnj^i^i " — ifi + 1) which can 
be absorbed into the normahzation of Z. The contribution for antifundamentals is 
obtained by the substitution mf — t- ei + e2 — mf. We note that the potentials for 
the fundamental matter include Penner-like factors of the form mf log(x — mf) which 
appeared in the work of Dijkgraaf and Vafa [21]. Such terms also arise in matrix 
models for 2* theories analyzed in 



3. As in \21\ |22| [28] . in the four- dimensional case we introduce the function 

e.(^) = -xT{-x)T{x)e— - ^""^'^^ 



sin(7rx) ' 

with simple poles at all integer values of its argument. Upon integration along the 
contour C encircling oo[ part of the real axis, this function can be used to pick 
up all integer values of h^P G oof. 

4. Finally we replace the summation over U by the integration over Xi according to (128|) . 
This leads to the following expression 

Z''^=l I)Me-i^^^"(^\ VM = \{dx,\{{x,-x,f^ (33) 

where M G MatnN represents nN x nN matrix form the /3-ensemble, and the measure 
VM involves the /3-deformed Vandermonde determinant in the leading order in 62. 
Similarly, to the leading order the potential is given by 

V'^ix) = tx + KtKx) + Vyi,(x; m) 

where we included the linear term arising from the instanton counting parameter A. 
The contributions from vector and fundamental multiplets are given respectively in 
(I3T]) and f l32|) . The example of a potential for pure SU{2) theory is given in figure [3l 



4.2 Example — SU{2) theory with 4 hypermultiplets 

SU{2) theory with 4 hypermultiplets is a simple and important example relevant for 
the AGT correspondence. Let us denote masses of two fundamental multiplets by mi, m2, 
and two antifundamental ones by —m^, —m^^. We also put ai = —02 = a. Then the matrix 
model for this theory is given by (133|) . with the potential which can be written as 

T/4d {x-af{x + af .w^^, 

^ * {X — mi){x — m2){x — m-i)[x — m^) x + a 
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Figure 3: Matrix model potential for the pure SU (2) Seiberg- Witten theory. 



4 

+ nif log(x — mf). 
f=i 

The logarithmic Penner-hke terms in the second hne appear similarly in [24j . As discussed 
in [23j, to obtain the instanton contribution to the SU{n) theory, apart from setting 
ai = —02 one still has to isolate the appropriate U{1) factor. It would be interesting to 
understand if such factors play any role from the matrix model perspective. Nonetheless, 
we suppose we can treat them as an overall contribution, without changing the form of the 
matrix model. 

5 Matrix models for five- dimensional theories 

In this section we derive matrix models for five- dimensional theories. To start with, we 
follow the steps 1—4 from the previous section, in case with mcs = 0. Next we discuss 
arbitrary mcs, which corresponds to turning on Chern-Simons terms. 

5.1 Theories with fundamental matter 

Here we again follow the steps 1 — 4 listed in the previous section, now in the context 
of five- dimensional theories. 

1. The factor (125]) for four-dimensional theory is realized in terms of f^'^{x) = {q^; q)^ 
functions (fT8|) . With analogous identification as in fl29|) . expressing the result in 



18 



terms of Xj variables fl26p . and using the asymptotics (157|) we get 

Up to the overall phase factor, this can be written as 

n(2sinh^i^) . (35) 



i<j 



2. The remaining factors from Z^^^, given in (12^ with f^'^{x) = {q^; q)^, can be writ- 
ten as in (p7|) . Using the notation of the quantum dilogarithm (15^ . and then the 
asymptotics (155]) . this leads to the expression 



n(n-h+bi-(imodN-N)(l-l3)\ n nN , . 

TT TT ^ ^) ^ rr TT p4 {Li,{e^^-^l)-U,(e-^.+'^l)) 

1111 0('0^>-bi + (WiV-Af)(l-/3)+n -1111 
«=1 i=l ^ 



e'2 



where Xi = €2/2- From the inversion relation for the dilogarithm 

_ 1 TT^ 

112(2;) +112(2; 1) = --(log2;)2 + _ -zvrlogz, 

we get, for the U{n) theory, up to constant and imaginary factors and to the leading 
order in €2 

n 

K'el(^) = ^^' + 2$^Li2(e-+"')- (36) 
1=1 

This potential agrees with results of [21jJ^I In a similar way, the contributions from 
fundamental multiplets fl25|) give 



^/w(^) = ELi2(e™). (37) 
f 

3. The quantum dilogarithm (15^ vanishes, giq'^) = 0, for h a positive integer. At such 
points its derivative is 

(1)2 iTvh -h{h-l)/2 

q'(q^) = -^^ ' ^ 

Therefore the following function has simple poles with residue 1 for x = q^ with 

hen 



fpolesi'^) 



(1 - x)v/a;5'(a;)5'(a;~^) 



^The opposite sign of x in the exponent appears due to a different convention q — e*^^ here, versus 



q = e 3= in 
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Written in the exponential form, it contributes only a linear term in x to the potential 

[211 Eg. 

4. Finally we get 

Now the measure is given by the Vandermonde determinant which is both /3-deformed 
and sinh-deformed, while the potential reads 

V'\u) = tx + V!i{x) + VfU^\ 

where we again included the linear term arising from instanton counting parameter 
A, as well as fpoiesi^)- The contributions for vector multiplet and hypermultiplet are 
given respectively in (j36|) and (137|) . 

5.2 Chern- Simons terms and more general Calabi-Yau manifolds 

In five- dimensional theories one can also include Chern-Simons terms \5^ ES] of the 
form 

j CijkAiAFj AFk, 

with indices labeling vector multiplets. It turns out that there is only a discrete consistent 
choice of Cijk and they are labeled by a single integer rrics- Five-dimensional theories with 
these Chern-Simons terms are also equivalent to topological string theory on appropri- 
ate toric Calabi-Yau manifolds [lOj. In this context the constants Qjk translate to triple 
intersection numbers of these manifolds. This dictionary was discussed in detail in [21j. 

The equivalence between Nekrasov partition functions, with ei = —£2 = h, and topo- 
logical string amplitudes computed from the topological vertex has been explicitly checked 
in [12]. It turns out that this equivalence extends also to the case of arbitrary ei, £2. It was 
proposed that topological string amplitudes in this case should be refined to two parame- 
ters [q = e'^^,t = e"*^^) in terms of the refined topological vertices [T^ [T5]. In |16[ [T7] it has 
been checked that such refined amplitudes indeed reproduce Nekrasov partition functions 
and the results for both vertices agree, and the consistent refinement of the Chern-Simons 
terms has also been found. It was shown that the refined Chern-Simons terms take form 
[161 [IT] 

1=1 
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where ||A|| = Yli ^i- This can be written exphcitly in terms of the lengths of rows 

1=1 

n 

1=1 

where n\ = Aj(Aj + 1 — 2i) = ||A|p — ||A*|p. Therefore, this expression can be written in 
terms of hf^ using ( IT9|) . 

Nonetheless, the resulting expression cannot be interpreted as a one-matrix model for 
non-zero mcs- This is so, because each set of hf^ is coupled to different ai through the 
term Therefore the symmetrization (j28|l cannot be performed for the entire set 

of Zj's, but only within each set /(fe-i)7v+i, • • • , IkN- In consequence we obtain the n-matrix 
model, with n sets of eigenvalues xf^ = e2hf\ and with the leading contribution to the 
potential of the linear form 

n N 

V''^' = mcsJ2{^^T.'^?)- (40) 

1=1 i=l 

We also note, that while it is straightforward to include the quadratic terms appearing in 
the last exponent in (1391) . they cannot be reinterpreted as the contribution to the potential 
of a matrix model due to an explicit dependence on i. Happily these terms are sub leading 
in 62- Therefore (HOj) is the only leading contribution to the matrix model potential V^^dcs^ 
and this is the same as in ei = — £2 = h case discussed in |21] . 



5.3 Relation to other deformations 

We found that in five dimensional theories, in the leading €2 expansion, the Vander- 
monde determinant takes a form of a sinh-deformation (135|) . familiar from |13[ [2T| |22| |28] . 
In the context of matrix models related to Nekrasov partition functions, another - seem- 
ingly unrelated - deformation was postulated in [47], which amounts to replacing the 
/9-deformed Vandermonde determinant by 

v=\{>^rii^T^- (41) 
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oo 



The identification A/ = g'^ makes contact with our notation. As discussed in [17] the 
integral^ 

S\a„a,,fi;z) = J fld,\iV,\z - Xjq-'/')Xf-'/^^ ^^^^ 

which includes this supposed deformed determinant, reproduces Nekrasov partition func- 
tion for SU{2) theory with four fundamentals, with a special choice of parameters 
0:3 = — ei/2. This expression is in fact a Jackson integral, i.e. a discrete integral de- 
fined by J dqxf{x) = (1 — g) 'YlT=o f i^^)^^ •> which reproduces Riemann integral in g — 1 
limit. The value of the above integral can be expressed in terms of g-deformed Jack or 
Jacobi polynomials. 

We now note several similarities between these expressions and matrix models discussed 
in this paper. Firstly, the discreteness in variables A/ in the Jackson integral is of similar 
kind as the discreteness of eigenvalues in our matrix models discussed in previous sections. 
Secondly, we observe that the deformation (HT]l is closely related to the determinant-like 
expression arising explicitly from rewriting Nekrasov partition function, given in 
Indeed, both expressions differ by terms of the form {imodN — jmodN){^ ~ P) in 
which from our perspective get multiplied by 62 and are subleading. Up these terms, the 
form of both expressions in the asymptotic form (157|) is the same, and in particular they 
both lead to the sinh-deformed Vandermonde (135|) in the leading order. The insertions 
of (A/g''^"^/'^^; g)oo in the integral 5'''(ai, 02, /3; 2;) are given by the same infinite products 
as contributions from vector- and fundamental multiplets in (12^ and (125|) with f^'^{x) = 
{x; q)^- This suggests, that one might indeed express matrix models for Nekrasov partition 
functions, to all orders in £2 deformation, in terms of Jack or Jacobi polynomials. 

Moreover, it was suggested in [17] that the deformation of the measure to the five- 
dimensional case could be related to the g-deformed Virasoro algebra considered in |44] . 
If this is indeed the case, it would be interesting to see how such an algebra manifests 
on the level of our matrix models. For example, in the standard formulation of matrix 
models, the loop equations are equivalent to Virasoro constraints. These constraints can 
be written in terms of operators which satisfy the Virasoro algebra and annihilate the 
matrix model partition function. It is interesting whether a similar structure arises in the 
g-deformed case. It is also tempting to extend such a deformation of Virasoro algebra to 
the six-dimensional case. 



4 



Note a different sign convention €12 — > £1,2 — —^1,2 in |47| . 
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Certain deformations of Virasoro algebra are also related to the so-called Q-bosons and 
Q-fermions, and Q-deformed boson-fermion correspondence. These objects also arise in 
the context of topological strings |56]. Even though in this case the role of the deformation 
is different, i.e. it plays a role of the Kahler parameter, it would be interesting to study if 
there are some relations between these both deformations of the Virasoro algebra. 

6 Further research 

In this paper we derived matrix models for /3-ensembles, which encode the instanton 
part of Nekrasov partition functions. This is just the first of several steps which should be 
completed. Most of all, it is important to analyze these models using matrix model tech- 
niques and confirm the relation to Seiberg-Witten theories from this perspective. Typical 
matrix model analysis involves finding the spectral curve and solving the loop equations. 
In the present case this is entirely non-trivial: as shown in [261 EZ], matrix models for 
/3-ensembles lead to non-commutative spectral curves. Moreover the form of our potentials 
is more complicated than the polynomial potentials analyzed in PSIEZ], and therefore the 
entire theory has still to be extended. Apart from general interest, the ultimate goal of 
this program, as explained in the introduction, would be to provide more rigorous proof of 
the AGT conjecture. 

Even though our models share several similarities (deformed measures, Penner-like 
potentials) with matrix models discussed in [211 HS| W7\ . there are also some differences. 
Apart from the Penner-like terms, our potentials include additional contributions. Some 
parameters arise on different footing: for example the instanton counting parameter A is 
encoded in the linear term in our models, while in [24j is appears under the logarithm. 
It would be important to elucidate these discrepancies. It is also important to extend 
our analysis to bifundamental fields and quiver gauge theories, and correspondingly Toda 
systems. Matrix models for such six- dimensional theories, with elliptic deformations, can 
also be written down explicitly. 

It would also be interesting to reexpress our results as matrix models in terms of Jack 
or Jacobi polynomials. Supposedly, the issue of symmetrization of the subleading terms 
should not arise in this case. Nonetheless, to the leading order such putative matrix models 
must agree with ours, and they should lead to the same (non-commutative) spectral curve 
and the entire solution. 

Assuming that the dualities discussed in the introduction indeed hold, it is important 
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to understand how various features of gauge theories and Liouville theories manifest in 
terms of our matrix models. One should understand the role of the one-loop part of the 
Nekrasov partition functions, as well as the decoupling of the U{1) factor, which play an 
important role in the AGT conjecture. Some other such features include: the relation to the 
g-deformed Virasoro algebra (discussed briefly in section [5^ . extension of this relation to 
the six- dimensional elliptic case (and finding appropriate more general deformation of the 
Virasoro algebra), interpretation of recursion relations between conformal blocks, matrix 
model interpretation of surface and loop operators discussed in |32l |l3], and many others. 

Finally, it is desirable to elucidate physical and mathematical interpretation of our 
matrix models. Physically, it is tempting to connect our results to phenomena involving 
D-branes and geometric transitions, underlying the whole Dijkgraaf-Vafa program [8|[9|l2^. 
Mathematically, our results hint into new structures hidden in instanton moduli spaces, 
which underlie the localization program of Nekrasov j31 IH El E]- 
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A More on Nekrasov partition functions 



In general, instanton parts of Nekrasov partition functions take the form 



^ ^ ^ ^ ^X,vec ^ \,{anti) fund ^ X,adj ^ X,bifund 



7CS 

A={A(i),...,A(")) 



where the summation extends over n sets of two-dimensional partitions A*-'-*, |A| = • ■* 
is the total number of boxes in a given set of such partitions, and A represents appro- 
priate instanton counting parameter (i.e. dynamically generated scale in asymptotically 
free theories, or appropriately renormalized bare coupling in conformal theories). Z^^^^, 
Zx,(anu)fund^ ^MP ^\,bifund ^nd Z?^^^ represent respectively contributions from vector 
multiplets, hypermultiplets in the (antifundamental or) fundamental representation, hyper- 
multiplets in the adjoint representation, in the bifundamental representation, and (trivial 
in four-dimensional theories) Chern-Simons terms. These terms depend on the Coulomb 
branch parameters ai and masses of matter fields rrif. The form of these terms depends 
also on the gauge group [7j. In this paper we focus on the U{n) gauge groups, so we 
write below explicit expressions only in this case. The form of all these terms depends also 
on the dimensionality of the spacetime: five- and six dimensional terms are respectively 
trigonometric and elliptic generalizations of the four-dimensional ones |5l|6l[T3]. The con- 
tribution for antifundamental hypermultiplets antifund related to that of fundamental 
ones -Z^X/und replacement of mf by ei + £2 — Trif, so usually we write down explicitly 
only the latter ones. The contribution from adjoint or bifundamental hypermultiplets is 
not considered to much extent in this paper. To find the total partition function one has 
to insert into the above sum the appropriate number of these terms, corresponding to the 
field contents of the theory of interest. 

There are several representations of Nekrasov partition functions which are somehow 
scattered through literature, see e.g. |3l HI O |6l [TTl [161 [23]. We discuss them briefly below. 
The representation which is most useful to construct our matrix models, discussed as the 
third one, appears to be used the least often. 

Firstly, various factors Z^ introduced above can be written as product over all boxes in 
partitions A*-'^. In case of vector, bifundamental or adjoint multiplets these products involve 
the arm-length ay^(i){\3) = Xf^ — j and the leg-length = X^^'^ — i (shown in figure [2]) 

of each box □ = G A*-'^ The explicit form of these expressions is nicely summarized 
e.g. in j23]. As they are not very useful for our purposes, we just just recall as an example 
that the contribution from fundamental hypermultiplets in the four-dimensional case takes 
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the form ^ ^ 

^xjw=n n («^+^ei+je2-m), 

1=1 □=(ij)GAW 

while vector multiplets in the five- dimensional case contribute [T7] 

n 

r^5d ^ TT ^ 

The second form of Nekrasov partition functions involves lengths of rows of partitions 
A*^'-*, as well as transposed partitions A'-'-*'*. In particular, the contribution for vector mul- 
tiplets reads [1] 



n 



ai - at + ei(z - 1) + e2(-j) 



while several equivalent forms of the five- dimensional contribution can be found in |17) . 
The contribution for Chern-Simons terms in five- dimensional theories (discussed in section 
15. 2p reads 

1=1 

and this is trivial in the four- dimensional case. 

The third form of Nekrasov partition functions involves just lengths of rows Xf^ and 
their differences. In such a form, the four-dimensional contributions for vector and funda- 
mental multiplets read [4J 

eT^'\JA.) r(Af -Af + /3(j-z) + M r(/3(j-0 + % + ^ 

^XJund 1111 Y{hi- M -i(3+l) ' ^ ' 

while those in five dimensions pTj Wl\ 

Zu d = TTTTTT^I (45) 

A,/W llUll^^,(0+,^_,,^_,^+,.^^^ 

The form of the five-dimensional terms is in fact subtle, and depends on the operator to 
which the localization formula is applied. The more often encountered form involving sinh 
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functions arises in case of the computation of the index of the Dirac operator, while the 
above form would correspond to the index of the Dolbeault operator. We find the above 
form more convenient, and the difference between the formulas involving sinh functions 
amounts simply to a redefinition of the linear term in the potential of matrix models which 
we find. From the geometric engineering viewpoint one could in fact consider various 
blow-ups representing matter multiplets, however we restrict only to the simplest case 
given above. 



B Rewriting the sums 

For an arbitrary function /, we wish to rewrite the contributions from the vector 
multiplet, given in (120]) 

„ ^ TT - fef + (' - + fjkk + I3U - i)) 

" <M,il., - + (' - m - /?)) + "(J -')+py ' 

in a form suitable for matrix model interpretation. To start with, we restrict to two sets 
of partitions, n = 2. In this case we denote the 2-tuple of partitions as 

(A,/i) = (A«,A(2)) = A, 

and instead of general (IT^ we introduce 

h.^h^'^ = X^-i + N + h, k, = hP = fLi-i + N + h. (47) 

Then (H^ takes the form 

Zn=2 = ZhhZkkZhkZkh- (48) 

where the consecutive terms arise respectively from the products in (HHj) over I = k = 1, 
I = k = 2, {I = 1, k = 2) and {I = 2, k = 1), which we consider now one by one. 
Firstly we consider Z^h, which reads 

and and correspond to products with, respectively, i < j and i > j. We consider 
first Z^f^ and split the product into three parts: 1 < i < j < N, {1 < i < N,N < j), and 
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N <i < j. This leads to 

f{h.-h, + {t-j){l-^) + P) 



Af oo 

X 



iiii /(/i,-6i + (i-7V)(l-/3)+j/3) + 

/J^ /(fe,-fe, + (i-i)(l-^)+/3) wJ^ ^(^) X 

\\}. f{K-hj + {i-j){l-P)) Jyi\f{h,-h + (i-N){l-P)+p))' 

where we used HJli ^/(t+j-^f^ = 7(i^- ^® 



7> _ f^ f{h^-h, + {^-J)il-P) + P) ^^^ f{-{h,-h + i^-N){l-(3))) ^ 
^f^f^ 1,11 . All f(o) ) 



i>j 



(50) 

The factor Zkk can be written in the same way as above, with hi and bi replaced 
respectively by ki and 62 



hi—^ki,bi—^b2 



Zkk = Z^k Zkk — Zhh Zhh 

Next we consider Z^ki which we write as 

The three factors Zj^jf^'^ on the right correspond respectively to i = j, i < j and i > j, 
and similar manipulations as above lead to the following results 

^< ^ / ^ / (/^. - % + (i - j) (1 - ^) + ^) w ^ + ^) \ 



+ All/(/^,-62 + (i-iV)(l-^)+/3) 

^ . /(/^, - fc, + - - /3) + /3) ^ ^ /( - (fc, - 61 + (z - iV)(i - 



Ul fr/,,-K + r,-7vi-/3ii All ff5^,) )■ 



f{h,-k, + {i-3){l-P)) nil /(5^^) 

Note that in the overall product of these three factors, all terms which involve /(612) and 
/(&12 + /5) cancel. 

The factor Zkh can be obtained from Z^k by the substitution hi k^ and 61 •<-> 62 



■^fc/i = -^fe/i -^4 -^fc/i - -^ftjt -^4 -^wt 
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hi-(r^ki,bl-(r^b2 



The expression fl48p is a product of the above factors, and from their form it is clear 
that it easily generalizes to arbitrary n: one has just to consider the analogous factors as 
above for all possible pairs (/, /c), with Z, /c = 1, . . . , n. This expression has also a simpler 
representation in terms of 

7 ._ ,'7,W 7,W 7,(1) 7,(1)^1 ^ 7 ._ 7,("-fc+l) 



In this notation, and denoting imodN = imodN, the expression (H6|) for arbitrary n takes 
the form 



Z 



AirT /(-(^.-fe^ + (w^-iV)(l-/3))) 
''fifi/(^.-&z + (W7v-iV)(l-/3)+/3)- ^ ^ 

C Asymptotics 

C.l Asymptotics of gamma function and quantum dilogarithm 

The following expansion of the logarithm of the gamma function holds |57] 

1 1 °° B 

logr(z) = -- log \z\ +zlog\z\-z + - log27r + 2n{2n -W--^ ' ^^^^ 

n=l ^ ' 

We also consider the asymptotics of the quantum dilogarithm, in the notation of |2I] 
written as 

00 ^ 

^W = n(i-i^O = (';^)oo' (54) 

where the g-Pochhammer symbol is 

oo 

For q = we have 

1 °° 1 R 

^ogg{z) = -J2u,_^{-)-^i-e,r- (55) 
€2 1 z m\ 

171=0 
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C.2 Asymptotics of the ratio of gamma functions 

The asymptotics 

+ (56) 



T{z + a + P) _ 



r{z + a) 

^ ' rt=l 

was derived in [58]. The first term of this expansion |57] can be obtained from the Stirhng 
formula, though to get the entire expansion requires more work. The coefficients are 
given by 



^-^ \n — m I 



1 """^ / 



m=0 

where coefficients Am satisfy the recursion formula 

n-l 

/ (T — m. \ 

^ \Arr,ia] 

\n — m, 4- 1 / 

m=0 

and we also set Cq = 1. From this prescription one can find 

One can also write 
where c„ are given in terms of the generating function 



= l/3(/3 + 2a-l), C^W^a) = ^ (2) ((/? - 2)(3/3 - 1) + 12a(/3 + a - 1)) 



etc. 



n=0 

C.3 Asymptotics of the ratio of quantum dilogarithms 

In the five- dimensional case we also need the following asymptotics 

— — (57) 

(^g"+^;g)oo 

which can be found explicitly e.g. in [59]. This formula is closely related to the ratio of 
g-gamma functions Tg{x). The function Tg is defined as 

^'<^>^<^->""^(?^- 

and it is known that 

lim r„(a;) = T{x). 

This function can also be expressed as the g-hyper geometric function. One can find the 
asymptotics (157|) using the g-analogue of the Stirling formula for g-gamma function, which 
is derived in |60| (see also |61]). 
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D Matrix model potentials 



As explained in the introduction, our methods allow to get the leading 62 terms of matrix 
model potentials. These leading contributions should be sufficient to get the spectral 
curve of the matrix model, and subsequently solve it. The subleading terms, for general 
/9, depend not only on xf\ but also explicitly on (imodN). Therefore they cannot be 
simply symmetrized and their matrix model interpretation is less clear. Nonetheless, for 
completeness, we write down the full form of the four-dimensional potential f l3T]) arising 
from Z^*^ . Up to the symmetrization of eigenvalues it reads: 



y4d 
' vec 



n 

» = ^^-2{x-ai + t2{i„^odN-N){l- 13)) + 



1=1 

+ {x-ai + e2{imodN - - /?)) x 

X [log {x-ai + e2{i„,odN - N){1 - 13)) + 

+ log {x-ai + e2(wiv - N){1 -(3) + 62/?)) + 
+^log {x-ai + e2{imodN - N){1 - /?)) + 
-— log {x-ai + e2{i„,odN - N){1 - 13) + €2(3) + 
-62(3 + €2/3 log {x-ai + t2{imodN -N){1- [3) + €2(3) + 

+ tr 2^(2^ - 1) I (a; - a, + e2{imodN -N){1- P)Y''' ^ 

1 



(5^ 



{x-ai + e2{imodN - N){1 - P) + €2/?)^* ^ 

This form arises from asymptotics presented in the appendix O) and contributions for other 
factors are analogous. 
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